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Abstract: Cubature formulas, asymptotically optimal with respect to accuracy, are derived for 
calculating multidimensional weakly singular integrals. They are used for developing a universal code 
for calculating capacitances of conductors of arbitrary shapes. 

1. Introduction 

Optimal with respect to accuracy methods for calculating singular integrals are being actively 
developed presently. They represent an important field of computational mathematics. Asymptotically 
optimal and optimal with respect to order ( to accuracy and to complexity) algorithms for calculating 
singular integrals on closed and open contours, and multidimensional singular integrals have been 
constructed in [1-3] on Holder and Sobolev classes of functions. 

In constructing optimal with respect to accuracy methods for calculating one-dimensional, bisin- 
gular and multidimensional singular integrals, a general method, proposed in monograph [1], was 
used. This method can be applied not only to singular integrals but also to other integrals with 
moving singularities. 
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This method allows one to construct several asymptotically optimal and optimal with respect to 
order and to accuracy algorithms for calculating hypersingular integrals [4], the Poisson and Cauchy 
type integrals [5], and multidimensional Cauchy type integrals. 

Although multidimensional weakly singular integrals are used in many applications, optimal meth- 
ods for calculating these integrals are not developed. 

An exception is the book [1], where asymptotically optimal with respect to accuracy methods for 
calculating integrals of the form 

2tt 2tt 




ct 9 ~ — 
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, CT2 - S 2 

ct 9 ~ 



72 



< 71, 72 < 1, were constructed on Holder and Sobolev classes. 

Thus, the development of optimal methods for calculating multidimensional weakly singular inte- 
grals is an actual problem. Construction of efficient cubature formulas for calculating weakly singular 
integrals for calculating capacitances of conductors of arbitrary shapes by iterative methods proposed 
in [6] and [7] is very important in many applications, for example, in wave scattering by small bodies 
of arbitrary shapes and in antenna theory. A bibliography on methods for calculating capacitances 
and polarizability tensors is contained in [7]. 

In this paper the method proposed in [1] is generalized to multidimensional weakly singular 
integrals. As a result the analogs of the basic results for singular integrals, obtained earlier, are 
obtained for weakly singular integrals. Moreover, we study the applications of optimal with respect to 
order cubature formulas for calculating weakly singular integrals on Lyapunov surfaces. Our results 
are used for constructing an universal code for calculating capacitances and polarizability tensors of 
bodies of arbitrary shapes. 

This paper consists of two parts. 

In the first part of the paper optimal methods for calculating integrals of the types: 

Kf=JJ f(^^ida 2 < Sl)S2<27r . (L1) 



and 



o o 
l l 



^" R^ -'^^ °<*<'- M 



-1 -1 

are constructed on Holder and Sobolev classes of functions. 



Our results for integrals (1.1) can be generalized to the integrals with other periodic kernels and 
functions. The development of cubature formulas for integrals (1.1) is of considerable interest because 
the results are applicable to integrals with weakly singular kernels defined on closed Lyapunov surfaces. 
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It will be clear from our arguments, that the results can be generalized to I— dimensional integrals, 
Z = 3,4,-- - . 

The second part of this paper deals with the iterative methods for calculating capacitances of 
conductors of arbitrary shapes. A general numerical method for calculating these capacitances is 
developed, and the results of numerical tests are given. 

2. Definitions of optimality. 

Various definitions of optimality of numerical methods and a detailed bibliography can be found 
in [8,9,10]. Let us recall the definitions of algorithms, optimal with respect to accuracy, for calculating 
weakly singular integrals. 

Consider the quadrature rule 

ni n 2 pi p2 

T /= E EEE?t.wii s (ti>#' ,8) fe,!fe)+ 

fc 1 =lfc 2 = l/ 1 =0«2=0 

+Rnin2(f;Pk 1 k 2 lih', x k 1 ,yk 2 ',tl,t2), (2.1) 

where coefficients p^fohh (*1j*2) and nodes (xfci>2/fc 2 ) are arbitrary. Here f( ll,l2 \si, S2) = d ll+l2 f (s\, S2) / ds 1 ^ ds 2 2 . 
The error of quadrature rule (2.1) is defined as 

Rn 1 n 2 (f',Pk 1 k2hl2' X k 1 ,yk 2 ) = SU P I Rn in2 (f', Pkik 2 hh ! x ki , Vk 2 ', h , h) \ ■ 

(tl,t 2 )6[-l,l] 2 

The error of quadrature rule (2.1) on the class ^ is defined as 

Rn 1 n 2 {^',Pk 1 k 2 hl 2 ^ x k 1 ,yk2) = SU P Rn in2 (/, Pkik 2 hl 2 ! ^fci , Vk 2 ) • 

Define the functional 

Cn in2 {^) = inf Rn 1 n 2 (^;Pk 1 k 2 i 1 i 2 ;x kl ,yk 2 )- 

Pk 1 k 2 l 1 l 2 i x k 1 ,Vk 2 

The quadrature rule with the coefficients Pk lk2 hl 2 anc ^ ^he n °des (a?^ , j/^ 2 ) is optimal, asymptot- 
ically optimal, optimal with respect to order on the class ^ among all quadrature rules of type (2.1) 
provided that: 

Rn 1 n 2 {'^'->Pk 1 k2hl2 lX *kl' i yW i 1^1 

The symbol a x j3 means Aa < (3 < Ba, where < A, B < 00. 
Consider the quadrature rule 



Tf = 5>fc(ti,* 2 )/(M fc ) 



fc=i 

+R n {f; Pk ;M k ;t 1 ,t 2 ), (2.2) 
where coefficients Pk{ti,t2) and nodes (M&) are arbitrary. 
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The error of quadrature rule (2.2) is defined as 

R n {f;p k ;M k ) = sup \R n (f;p k ; M k ;t 1 ,t 2 )\. 

(tl,t 2 )6[-l,l] 2 

The error of quadrature rule (2.2) on the class ^ is defined as 

R n (V;p k ;M k ) = sup R n {f,p k ; M k ). 



Define the functional 



CnW = mf Rn(^; Pk ;M k ). 

Pk\M k 



The quadrature rule with the coefficients p* k and the nodes (M£) is optimal, asymptotically 
optimal, optimal with respect to order on the class \I> among all quadrature rules of type (2.2) provided 
that: 

R n (*;pt;M*) 

cm = 1 '~ 1 '~ 1 ' n ^°°- 

By R nin2 (^) the error of optimal cubature formulas on the class * is defined. It is obvious that 

Rn in2 W = Cmn 2 (*)- 

3. Classes of functions 

In this section, we list several classes of functions which are used below. Some definitions are 
from [11,12]. 

A function / is defined on A=[a,b] or on A=K, where K is a unit circle, satisfies the Holder 
condition with constant M and exponent a, or belongs to the class H a (M), M > 0, < a < 1, if 
\f(x') - f(x")\ < M\x' - x"\ a for any x',x" G A. 

Class Hu, where u>(h) is a modulus of continuity, consists of all functions / G C(A) with the 
property |/(xi) - f(x 2 )\ < Mu(\xi - x 2 \),x 1 ,x 2 G A. 

Class W r (M) consists of functions / G C(A) which have continuous derivatives /', /", . . . , /( r_1 ) 
on A, and a piecewise-continuous derivative on A satisfying max xe [ a ^\f^ r \x)\ < M. 

Class Wp(M), r = l,2...,l<p<oo, consists of functions f(t), defined on a segment [a,b] or 
on A = K, that have continuous derivatives /', /", . . . , f^'~ l \ and an integrable derivative such 
that 

i/p 

\f( r \x)\ p dx < M. 

lA 

Class Wa(M), r = 1,2 . . . ,0 < a < 1, consists of functions f(t), defined on a segment [a,b] or on 
A = K, which have continuous derivatives /', /", . . . , f( r \ such that 

|/M(zi)-/< r )(z 2 )| <M\ Xl -x 2 \ a . 

A function f(x±,x 2 , . . . , xi), I = 2,3, ... , defined on A = [tti, b 2 ;a 2 ,b 2 ; 
. . . ; ai, bi\ or on A = K\ x K 2 x ■ ■ ■ x K\, where Ki, = 1, 2, ... ,1, are unit circles, satisfying Holder 



conditions with constant M and exponent CKj, i = 1, 2, ■ ■ ■ ,1, or belongs to the class -ff Q , li ... jQi (M), M > 
0,0 < a < l,i = 1,2,..., I, if 

|/(xi,x 2 , ■ ■ • ,a;j) - /(yi,y 2 , • • • ,yi)\ < M(|xi - Vl \ ai + ■ ■ ■ + \x t - Vl \ ai ). 
Let to, Ui, where i = 1, 2, . . . , I, I = 1, 2, . . . , be moduli of continuity. 

Class H UljmmmjWl (M), consists of all functions / € C(A),A = [a\, b 2 ; a 2 , b 2 ; • • • ; a/, 6/] or A = K\ x 
K2 x ■ ■ ■ x Ki, with the property 

|/(xi,x 2 ,.. • ,xj) - f(y 1 ,y 2 ,.. ■ ,2/01 < M(wi(|xi — yi|) H hu;;(|x; - y ; |)). 

Let Hf(A),j = 1,2,3, A = [ai, 6 2 ; 02, ^2; • • • ; aj, 6j] or A = i^i x if 2 x ■ ■ ■ x K u I = 2, 3, . . . , be the 
class of functions f(x\, X2, ■ ■ ■ , xi) defined on A and such that 

\f(x)-f(y)\ <co( Pj (x,y)),j = 1,2,3, 

where x = (xi, . . . , x t ), y = {yi, . . . ,yi), pi(x,y) = maxi<i<i(\xi - yi\), p 2 {x,y) = Y!i=i N ~ 

Vi\,pz{x,y) = ELi \ x i - Vi\ 2 } 1/2 - 

Let Hf(A),j = 1,2,S,A = [a 1 ,b 2 ;a 2 ,b 2 ; . . . ; at, k] or A = K x x K 2 x ■ ■ ■ x K h I = 2, 3, . . . , be the 
class of functions f(xi,x 2 , . . . , xi) defined on A and such that 

\f(x)-f(y)\<( Pj (x,y)r,j = 1,2,3. 

More general is the class Hpj(A),j = 1,2,3. It consists of all functions /(x) which can be rep- 
resented as /(x) = p(x)g(x), where g(x) € H?(A),j = 1,2,3, and p{x) is a nonnegative weight 
function. 

Let Z^{A),j = 1, 2, 3, be the class of functions f(x±, x 2 , . . . , x/) defined on A and satisfying 
|/(x) + f(y) - 2/((x + y)/2)\ < u{pj(x, y)/2),j = 1, 2, 3. 

Let Z?(A),j = 1, 2, 3, be the class of functions f(x±,x 2 , . . . ,xj) defined on ^4 and satisfying 
|/(x) + - 2/((x + y)/2)| < (pj(x, y)/2) a ,j = 1,2,3. 

Class Zpj(A),j = 1,2,3, consists of all functions /(x) which can be represented as /(x) = 
p{x)g{x), where g(x) G Z?(A),j = 1,2,3, and p(x) is a nonnegative weight function. 

Let W ri '-' ri (M),l = 1, 2, . . . , be the class of functions /(xi, X2, . . . , x/) defined on a domain A, 
which have continuous partial derivatives 

d\ v \f(x 1 ,...,x l )/dx r £ ...dxf ,0 < |v| < r - 1>| = Vl + ••• + vi,n > v t > Q,i = 1,2,..., I, r = n + 
■ ■ - + ri and all piece-continuous derivatives of order r, satisfying ||3 r /(xi, . . . , Xi)/dx r ^ . . . dx\ l \\c < M 
and \\d r *f(0, . . . , 0, x t , 0, . . . , 0)/<9x[ ! || c < M, i = 1, . . . , I. 
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Let Wp 1 ''"' 1 " 1 (M), I = 1, 2, . . . , 1 < p < oo be the class of functions f(x\,x 2 , . . . ,xi), denned on a 
domain A = [ai, &i; . . . ; ai, b{\, with continuous partial derivatives 

<9H/( Xl) ... ^ti/dxl 1 ...dx^',0 < \v\ < r - 1, \v\ = vi + ■ ■ ■ + v h n > vi > 0,i = 1,2, . . . ,l,r = 
r± H + ri , and all derivatives of order r, satisfying 

\\d r f( Xl , Xl )/dxYdxZ . . . dx? \\ LM) < M, 

Wl+V2+ ... +Vlf{xiA Q)/dx r ldx v 2 dx v t \\ Lp{[aiM]) < M , \v 2 \ + H + ■ ■ ■ + h| < r - n - 1; 



ll^ 1+ - + ^- 1+r 7(o, . . . , o, Xl )/dx^dx% . . . dxl^dx? \\ Lp([aiM) < m, \ Vl \+\v 2 \+- ■ < r-n^-i. 

Let A = [ai, b 2 ; a 2 , b 2 ; ■ ■ ■ ; a;, b\\ or A = K\ x K 2 x • • • x K|. Let C r (M) be the class of functions 
f(xi, x 2 , . . . , xi) which are defined in A and which have continuous partial derivatives of order r. 
Partial derivatives of order r satisfy the conditions 

..dWf{ Xl ,...,x l ) l 



v t IIC 



< M 



11 dxi 1 ...dx 1 ; 1 

for any v = (vi, . . . , w/), where Vi > 0, i = 1, 2, . . . , / are integer and ^!=i v i = r - 
By ^ we denote the set of periodic functions of the class 

It is known [13] that Lyapunov spheres are defined as regions bounded by a finite number of 
closed surfaces satisfying the three Lyapunov conditions: 

1. At each point of the surface a tangent plane (and, therefore, a normal) exist. 

2. If 6 is the angle between the normals at the points mi and m 2 , and r is the distance between 
these points, then 

6 < Ar x , < A < 1, 

where A and A are positive numbers which do not depend on mi and m 2 . 

3. For all points of the surface, a number d > exists such that there is exactly one point at 
which a straight line, parallel to the normal at the surface point m, intersects the surface inside a 
sphere of radius d centered at m. 

Let S be a Lyapunov sphere, and iV be the exterior normal to this sphere. We introduce a local 
system of Cartesian coordinates (x,V,(), whose origin is located at an arbitrary point too of S, the 
( axis is directed along the normal iVo at the point too, and the x an d V axes lie in the tangential 
plane. In a sufficiently small neighborhood of too, the equation of the surface S in the local coordinates 
(Xj V> C) nas the form 

( = F( X ,V)- 

Definition 4.1. [13] The surface S belongs to the class Lk(B,a) if F(x,r]) G W£(B), and the 
constants B and a do not depend on the choice of the point toq. 



4. Auxiliary statements. 

We need the following known facts from the theory of quadrature and cubature formulas. These 
facts can be found, for example, in [11], [12], [14], [15]. 

Lemma 4.1. Let *i be the class of functions W£(l), 1, 2, ... , 1 < p < oo, < t < 1, f(t) G 
and the quadrature rule 

\ n 

/ f(t)dt = Y,Pkf(tk) + Rn(f) 
fc=1 

be exact on all the polynomials of order up to p— 1, and has error i? n (\I/i) on the class ^i. Let ^ be 
the class of functions W£(l), r = 1, 2, . . . , 1 < p < oo, a < x < b, and g(x) G VFp (1). T/ien £/ie 
quadrature rule 

h r n 

/ = (b - a) ^2p k g(a + {b - a)t k ) + R n {g) 

a k=l 

has error R n {^2) on the class of functions ^2 and 

Rn{^2) = {b-aY +1 - l ^R n (^ 1 ). 



Theorem 4.1. [11] Among quadrature formulas 



J f( X )dx = Y,Y.puf {l) ^ + R (f) = L ^ + R(f) 



k=l 1=0 



the best formula for the class W£ (1) (1 < p < 00) with p = r — 1 and r = 1, 2, • • • , or p = r — 2 and 
r = 2, 4, 6, • • • , is the unique formula defined by the following nodes x* k and coefficients p* kl : 



,* k = h{2{k-l) + [R rq {l)] 1 ' r ), k = 1,2,... 



m, 



2h 2v+1 ( 
(Z = 0,l,...,p), p* kt2v = — ^"^(l), fc = 2,3,...,m-l; v = 0,l,... 



r — 1 



P*k,2v+i = I A; = 2,3, . . . ,m - 1; t> = 0, 1, . . . , 



r - 2 



, h = 2 _1 (m - 1 + [R rq {l[ 



ii/H-i 



r-l 



and R rq (t) is the Chebyshev polynomial t r + ^ deviating least from zero in the norm L q (—1, 1), 
where p^ 1 + q~ l = I. Here 



uw;(i)} = R n [w;{i) 



Rrq(l) 



2 r r\fFq~TT(m - 1 + [12 r ,(l)] 1 A 



Let a function /(x, y) be given on a rectangle D = [a, b; c, d]. Consider the cubature formula 

r. „ m n 

J J f(x,y)dxdy = ^2^2p ki f(x k ,yi) + R mn (f), (4.1) 



fc=i j=i 



defined by a vector (X, Y, P) of a nodes a < x\ < x 2 < • • • < x m < b, c < y\ < y 2 < • • • < y n < d, and 
coefficients pki- 

Theorem 4.2 [11]. Among all quadrature formulas of the form of (4-1) the formula 

,. ,. m n 

/ / f(x, y)dxdy = 4hqY,Yl f( a + ( 2fc " c + ( 2i ~ !)</) + ^(/), 
J J fe=l i=l 

where h = q = is optimal on the classes H^^^D) and H^{D). In addition 

h q 

Rmn[H^ 2 {D)\ =4mn[q J ^(fydt + h j uj 2 (t)dt}; 


q h 

R mn [H%(D)} =4mn J J uj^t 2 + T 2 )dtdr. 
o o 

Consider the cubature formulas of the form: 

JV 

J J p(x,y)f(x,y)dxdy = J> fc /(M fc ) + R(f), (4.2) 

D k=1 

where p(x,y) is a nonnegative and bounded on D function, p k , M k (M k G D) are coefficients and 
nodes. 

Theorem 4.3 [11]. Let p(x,y) be a nonnegative bounded weight function. If R]y[Hp j(D)] and 
R^[Zpj(D)], where j = 1,2,3, and < a < 1, are the errors of optimal formulas as (4-2) on the 
classes - {D) and Z'* j(D), respectively, then 



lim N a > 2 R N [H? AD)} = lim N a / 2 2R N [Z? AD)] = 

N^oo P ' J N^oo ' ' ' ■' 

- (2+a)/a 

j j(p(x,y)) 2 /^dxdy 



, j = 1,2,3, 



/ \ (2+a)/a 7r / 6 

where D l = ^ (^) / D 2 = 2 1 "«/(2 + «), and D 3 = 2 1 -/ 2 /(2 + a). 

7/j = 2, then the conclusion holds for n— dimensional cubature formulas. 
Remark. Theorem 4-2 is generalized to the case of unbounded weights p(x , y) in [2]. 
In this paper we will use the following result (see [16]): 

Lemma 4.4 . Let H be a linear metric space, F be a bounded, closed, convex, centrally sym- 
metric set with center of symmetry 9 at the origin, and £(/), l± (/),..., In(/), be some linear func- 
tional. Let S(li(f), . . . , Zjv(/)) be some method for calculating the functional L(f) using functionals 
(h(f), • • • , £jv(/)), and S be the set of all such methods. Then the numbers Di, . . . , D^ exist such that 

N 

sup \L(f) - V D k l k (f)\ = inf sup \L(f) - S(h(f), l N (f))\. (4.3) 
feF ^ s feF 



This means that among the best methods for calculating functional L(f): 

L(f)^S(h(f),...,l N (f)), (4.4) 

there is a linear method. 

Proof. Let us associate with each / G F a point (L(f), h(f), ■ ■ ■ ,In(I))- Let Y be a set of all 
such points (y , . . . , Vn) for / € F. 

From our assumptions, it follows that Y is a closed centrally symmetric set with the center of 
symmetry at the origin. 

Let (uq, 0, . . . , 0) be an extremal point of the set Y, and 

D = sup z = y . 

(z,o,-,o)eY 

Because F is bounded, one has Do < oo, and because F is convex and centrally symmetric with 
respect to the origin, one has Do > 0. 

Draw the support plane for the set Y through the point (Do, 0, . . . , 0) : 

N 

(yo-Do) + J2Cm = 0. 

3=1 

Since Y is centrally symmetric with respect to the origin, the plane 

N 

(yo + D ) + J2 C jyj = 

i=i 

is also a support plane for Y, and Y lies between these two planes. 
Hence, we have for the points of Y the inequality: 

N 

\yo - ^2 D jVi\ < D , Dj = -Cj. 

3=1 



The definition of yi implies 

N 



sup|L(/)- J>^(/)|<A). (4.5) 

Let /o be an element F corresponding the point (Do, 0, . . . , 0). Then S(h(±fo), ■ ■ ■ , Zjv(±/o)) = 
5(0, . . . , 0). The right-hand side of (4.3) is not less than 

inf max |L(/ ) - 5(0, ....,0)| = inf max{|Z? - a\, \D + a\} = D . 

S a 

This and (4.5) imply that the right-hand side in (4.3) is not less that the left-hand one. But the 
right-hand side of (4.3) can not be more than the left-hand side of (4.3) because a set of methods 5 
contains linear methods. Lemma 4.4 is proved. ■ 
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Corollary. Among all functions for which the optimal method for calculating L(t) has the greatest 
error for a given set of Junctionals, there exists a function satisfying the conditions h(f) = ••• = 
In(I) = 0. It follows from the proof that such a function is the function /q. 

5. Optimal methods for calculating integrals of the form (1.1). 

5.1. Lower bounds for the functionals Q nm and (jv- 

In this Section we derive lower bounds for the functionals Cnm and (n, defined in Section 2 , for 
calculating integrals (1.1) by the cubature formulas 

ni n 2 Pi P2 
fc 1= l fc 2 =l h=0 h=o 

+Rn 1 n 2 (f;Pk 1 k 2 l 1 h;X kl ,X k2 ;Si,S 2 ), (5.1) 

and 

N 

Kf = J2pk(s!, s 2 )f{M k ) + R N (f; Pk ; M k ; 8l , s 2 ) (5.2) 



k=i 

on Holder and Sobolev classes. 

Theorem 5.1. Let ^ = H UljU2 (D) or ^ = H^(D), and calculate integral (1.1) by formula (5.1) 
with pi = p 2 = 0. Then the inequality 

h q 



Cmrwt*] > J U>i(t)dt + h J L0 2 {t)dt], 



where q = — , h = — , and 

2tt 2tt 

f f ds\ds 2 t 

7 := J J {sin^{ Sl /2) + sin?{s 2 /2)Y ( } 



is valid. 

Corollary. Let <I> = H aa (D) or ^> = H^(D), and calculate integral (1.1) by formula (5.1) with 
n\ = n 2 = n and pi = p 2 = 0. Then the inequality 

, rTl 2 l^ a 

(1 + a)n a 

is valid. 

Proof of Theorem 5.1. Denote by ip(si,s 2 ) a nonnegative function belonging to the class 
#wiw 2 (l) an d vanishing at the nodes (x kl ,x k2 ), 1 < k\ < n±, 1 < k 2 < n 2 . 
One has: 

Rn\n 2 (V'i Vk\k 2 ! X kl , X k2 ) > 

2tt 2tt / 2tt 2tt 

> t~9 / / 1 /7r-^ — H 2: a2)da \ da2 — ^\d SldS2 = 



47r 2 i J \J J [sin 2 ((a 1 -s 1 )/2) + sin 2 ((a 2 -s 2 )/2))y 

\0 
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2tt 2tt / 2?r 2tt 

±J J^a 2 ) Iff [mn2{(ai _ Si)/2) dS + ld 2 2((f72 _ S2)/2 ))]A I ^ 
\o 

2tt 2tt 2tt 2tt 

Iff ds\ds 2 



j j tp(s 1 ,s 2 )dsids 2 . (5.3) 



4tt 2 7 7 [sin 2 (si/2) + sin 2 {s 2 /2)] x 



From Lemma 4.4 and Theorem 4.2 one concludes that the following inequality 

h 



7 

7T Z 



g y uji(t)dt + h j u 2 (t)dt 



7T 7T 

,h = —, g = — 
ni n 2 



o o 
holds for arbitrary weights Pfcife 2 an ^ nodes (x^^x^) and 

h q 

[g Jui(t)dt + h J ui 2 {t)dt\. 
o o 

Theorem 5.1 is proved. ■ 

Theorem 5.2. Zei ^ = Hf or ^> = Zf, i = 1,2,3, and calculate the integral Kf by cubature 
formula (5.2). Then 

( N [H? } = 2( N [Z?] = (1 + o(l)) 7 (47r 2 ) 2 / a AiV- a/2 , 
Proof. The proof of Theorem 5.2 is similar to the proof of Theorem 5.1, with some difference is 

2?r 2vr 

in the estimation of the integral f f ip(si, s 2 )ds±ds 2 , where the function ip(si, s 2 ) belongs to the class 

o o 

Hf (or Zf\ is nonnegative in the domain D = [0, 27r] 2 , and vanishes at N nodes k = 1, 2, . . . , N. 
Using Lemma 4.4 and Theorem 4.3, one checks that the inequalities 

2tt 2tt 



inf sup f [ i/}(s 1 ,s 2 )ds 1 ds 2 = (l + o(l))Di(4ir 2 ) {2+a)/a N- a/2 




2w 2tt 



inf sup f [ ij( Sl ,s 2 )d Sl ds 2 = (l + o(l))-Di(47r 2 )^ +a ^ a N' a / 2 

hold for arbitrary M k G D, k = 1, 2, . . . , N. 

Substituting these values into inequality (5.3), we complete the proof of Theorem 5.2. ■ 
Theorem 5.3. Let ^ = C 2 (l), and calculate the integral Kf by formula (5.1) with p\ = p 2 = 0, 

and n\ = n 2 = n. Then 

U[*]>(l + o(l))^, 

where K r is the Favard constant. 
Proof. Let 

i ! (si,s 2 ) = V'i(si) + ^2(52), 
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where < ipi(s) G VF r (l) vanishes at the nodes Xk, k = 1,2, ... ,n, and < ip2(s) £ W r (l) vanishes 
at the nodes yt, k = 1, 2, . . . , n. 

According to [11], for arbitrary nodes Xk, k = 1, 2, . . . , n one has: 

2tt 



/ 



i>i{s)ds > C —^ JL i i = l,2. 



Thus, the inequality 

> , , 8w 2 K r 

tp{Si, S2)dSidS2 > 



if 





holds for arbitrary nodes (x±, . . . , x n ) and (y\, . . . , y n ). 

The conclusion of Theorem 5.3 follows from this inequality and from (5.3). ■ 

Theorem 5.4. Let ^ = Wp T (l), r = 1, 2, . . . , 1 < p < oo, and calculate the integral Kf by 

formula (5.1) with p\ = p2 = r — 1 and n\ = ri2 = n. Then 

2 l l q ir r ~ 1 / p R (Y) 
C "" [ * 1 S (1 + ° (1)> r!(r, + !)■/,(„ 

where R rq {t) is a polynomial of degree r, least deviating from zero in L q ([—1, 1]). 

Proof. Let L = [j^]. Take an additional set of nodes £fc = 2 j^, k, I = 0, 1, . . . , L — 1. 

By (vi,Wj), i,j = 0,1, ... ,N — I, N = n + L, denote the union of the sets (xk,yi) and (£«,£?)• 
Let V'(si,S2) = Y>i(si) + Y>2(s2), where ^iX 5 ) £ Wp(l) vanishes with its derivatives up to the order 
r — 1 at the nodes Vi, i = 0, 1, . . . , N — 1, and ip2{s) S Wp^(l) vanishes with its derivatives up to 
order r — 1 at the nodes to,-, j = 0, 1, . . . , N — 1. Assume J ipi(s)ds > 0, i = 0, 1, . . . , iV — 1, and 

/ ip 2 (s)ds > 0, j = 0, 1, . . . , iV — 1, where ujv = 2ir and = 2ir. 

Wj 

Let 

0, i/ (0-1,0-2) = 0i,s 2 ), 



h(si,s 2 , 0-1,0-2) := 



! , otherwise, 



(sm 2 ((cri-si)/2)+sm 2 ((o-2-s 2 )/2)) A 



^ + (si,s 2 ) = < 



ip(si,S2), if ip(s 1 ,s 2 )>0, 
0, i/ ^(si,s 2 )<0, 



V 7 (si,s 2 ) = < 



0, */ ^(si,s 2 )>0, 
-ip(si,s 2 ), if ip(s 1 ,s 2 ) < 0. 
For each value i, j = 0, 1, . . . , N — 1, we have (with N = N\ = N 2 = L): 

2tt 2tt 



y y h(^ i ,^ j ,a 1 ,a2)ip(o- 1 ,a2)da 1 da 2 





AT-1AT-1 



^2^2 / ^(^' 0' cri, cr 2 )-0(cri, cr 2 )(icri<icr2 

J n 7 n ^ J 
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JV-1 JV-1 

^2^2 / ^(^'^' cr i' cr 2)V' + (ci ) cr 2 )dcri dc7 2 - 



fc=0 J=0 t , 
JV-1 JV-1 ^ fc+1 



^2^2 J J H&i^j^i,^)^ {^i,<J2)duida 2 > 



k=0 1=0 ? fe e, 



i+[(iVi-l)/2] j+[(JV 2 -l)/2] £*+i ft+i 

> X X] fr(&>€j>£fc+i>&+i) / / ip + (cri,cr 2 )daida 2 + 

k=i+i i=i+i £ ^ 

i+[(ATl-l)/2] j-1 & 

+ X] X] fr(&»£j,£fc+i,6-i) / / ip + (crua 2 )da 1 da 2 + 

k=i+l I=i-[(JVa-l)/2] £ ^ 

i-i j+[(jv 2 -i)/2] e* ft+i 

+ X] X] MC»)C?,6-l,Ci+l) / / V' + (o-l,CT2)'icrido-2 + 

fc=»-[(jvi-i)/2] j=j+i £ 

j-1 j-1 ^ & 

+ X] X] ^(&>£.?,£fc-l,£z-l) / / ^ + (0"l,0-2)^i(icr2- 

fe=i-[(JV 1 -l)/2] /=j-[(JV 2 -l)/2] 7 7 



i+[(JV!-l)/2] j+[(JV 2 -l)/2] &+i 



k=i+l l= i+ l ffc 5; 



i+[(JVi-l)/2] j-1 ^ 6 

XI X M&>€?>£fc>&) / / i/}~(<T 1 ,a 2 )da 1 d(72- 

k=i+l I=j-[(JV a -l)/2] £ ^ 

i-i j+[(jv 2 -i)/2] & 6+i 

X X ^(&> / / V 7 " (0-1,0-2) daida 2 - 

k=i-[( Nl -i)/2] 1=3+1 £ 

j-i j-i &f ^ 

X X ^(&> / / ip~((7i,a 2 )d(7 1 da2 = 

fc=i-[(iVi-l)/2] *=j-[(JV 2 -l)/2] f /_ i ^ 

i+[(JVi-l)/2] j+[(JV 2 -l)/2] €fc+l 6+1 

= ^ X ^te^i^fc+1,6 + 1) / / 1p((Tl,<T2)d(T 1 d<T2 + 

k^+1 1=3+1 I I 

i+[(JV!-l)/2] j-1 &+1 & 

+ ^ X ^(^,^,6+1,6-1) / / l/)((Tl,<T2)d(7ld<72+ 

k=i+l l=j-[(N 2 -l)/2] ( k ^ 

i-i j+[(JV 2 -i)/2] & ft+i 

+ X X] M&>£j,&-i,6+i) / / ^(o-i,a 2 )do-idcT2+ 

fc=i-[(ATi-l)/2] i=j+l £ 

j-1 j-1 ^ 

+ X X fr(&>€j>£fc-i,6-i) / / ^(o-i,(T 2 )dcri(io-2- 

fc=i-[(JVi-l)/2] *=j-[(JV 2 -l)/2] x ^ 
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i+[(iVi-l)/2] j+[(JV 2 -l)/2] 

Yl Yl (M6>^&>C0-M&>&>&+i>6+i))>< 

k=i+l l=j+l 

6+1 6+1 

6 6 

i+[(JV!-l)/2] j-1 

Y Y (M6,£j>6b6) ~ M6,£j>6fc+l,6-l))> 

k=i+l l=j-[(N 2 -l)/2] 

6 6 
6+1 6-1 

i-1 j+[(iV 2 -l)/2] 

(M&>&>&>£0-M6>&>&-i»6+i)) ; 

fe=i-[(ATi-l)/2] Z=j+1 

6 6+1 

x 

Sfc-i 6 

i-l i=j— 1 

X] (fr(6&,6,6)-M6>&>6fc-i,6-i)): 

fc=i-[(JVi-l)/2] J=j-[(iV 2 -l)/2] 

6 6 

x y y ip~(ai, cr 2 )daida 2 = 

6-1 6-1 

= Ji + J 2 + J 3 + Ja + h + h + h + h- 



y y ip~(cr 1 ,a 2 )da 1 da2- 



Let us estimate the integral 

6+1 6+1 



6+i 6+i 



y y V (cr 1 ,a 2 )da 1 da 2 < j J \ip (a 1 ,a 2 )\da 1 da 2 < 



6 6 



6 6 
6+i 



< 



6+i 6+i 6+i 6+i 

y y i^(<7i,<72)i<frid*2 < (6+1-6) y 1^)1^+ (6+1-6) y \m^)\^< 



6 6 



< 2 



lJ r!' 



where we have used the fact that the functions tpi(s) and if) 2 (s) on the segments [6)6+1] an d [6)6+1] 
vanish with derivatives up to order r — 1. 
Now let us estimate the sum: 

i+[(L-l)/2]j+[(L-l)/2] 

Y Y IM6, £?,6+i,6+i) - ^(6,0,6,6)1 = 

fe=i+i i=i+i 



15 



i+[(L-l)/2]j+[(L-l)/2] 

E E 



k=i+l 



=3+1 



stn z — 



(^) +sin 2l0±lzl)^ 



< 



< 



L 



< 



i+[(L-l)/2] j+[(L-l)/2] 

E E 



fe=i+l /=i+l 
i+[(L-l)/2]j+[(L-l)/2] 

; E E 

fe=i+i i=j+i 



(■ ovr(fc— i) . ■ o7r(/— i)\ 
sin — - + sm z V L y I 



l+A 



A; — i Z — 7 
+ 



L L 



< 



^2+2A 



( (A; _ i) 2 + (/ _ i) 2 ) l+A 



L 



< 



,'[(L-l)/2] [(i-l)/2] . 

< C (^) 2A E i+ E 



< c(L) 2A <( 



1 k=l 

logL if A = i, 
1 i/ A>± 



By c > various estimation constants are denoted. Thus 

1 



h=o 



n' 



The expressions 12,^3, and I4 are esimated similarly. 

From the definition of the function Y>(si,S2) it follows that the error of cubature formula 
for si = £j, S2 = £j can be estimated as follows: 

2w 2w 

R(ip,,ti,tij) = J J ^(vi,<T2)h(£i,£j,a 1 ,a 2 )daida2>o(J^J + 




i+[(L-l)/2]j+[(L-l)/2] 

+ E ^(&>£.?>6;+i>£h-i) / / i/}((Ti,a2)d(7 1 da2+ 

k=i+i l=j+i : ' ' 



?fe+i & 



i+[(L-l)/2] j-1 

+ ^2 E ^(&>£.?,£fc+l,6-l) / / 1p(<Tl,(T2)daid(72 + 

k= l+ l l=j-[(L-l)/2] I ^ 

i-1 j+[(L-l)/2] & 

+ X] E fr(&>€?>£fc-l,6+l) / / 1p(<7l,(72)d<Tld(T2 + 

k=i-[(L-l)/2] 1=3+1 I 

i-1 3-1 % % 

+ ^2 E ^(&>£.j>£fc-i,&-i) / / i\){a\,a-i)do\doi- 

k=i-[(L-l)/2] l=j-[(L-l)/2] ^ ^ 
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Averaging the above inequality over i and j, one gets: 



Rnn[$] > SUp max Rnnilp, > 



j L-li-l 

^EE 



i=0 j=0 



+[(L-l)/2]j+[(L-l)/2] 

E E 

fe=i+i i=j+i 
i-i 



5fc+i 6+1 



E E H€i:tj:€k+i,til+i) J J ^{cri,a 2 )duida2+ 



6 6 



i+[(L-l)/2] 

+ E 

fc=i+l i=j-[(L-l)/2] 



5fc+i 6 



E fr(&»£?,£fc+l>&-l) y y V ; (o"l,0-2)rfo-i(icr2 + 



€fc 6- 
6c 6+1 



i-1 j+[(L-l)/2] 

+ E E fr(&>£?>£fc-i>£/+i) / / ^{<ri,a 2 )daid(72+ 

fc=i-[(L-l)/2] /=j+l 



+ E E / / ip((Ti,(T2)d(7 1 da2 +o(— J 

fc=i-[(L-l)/2] i=j-[(L-l)/2] V " ' 



6-i 6 
6 6 



6-i 6-i 



>o - + 



n' 



1 



L-l L-l 



i+[(L-l)/2]j+[(L-l)/2] 

E E / / i>{?i,<72)<foid02 E E /i ^'^''^+ i '^+ i ) + 



fc=i+l l=j+l 



6 6 



i=o i=o 



L-l L-l 



i+[(L-l)/2] i=j-i 6+1 6 
fc=i+l i-[(L-l)/2] / fc ^ 



i=0 j=0 



L-l L-l 



i-1 j+[(L-l)/2] 6 ^+1 ^ , . , 



+ E 

fe=i-[(L-i)/2] i=j+i £ 



i-1 



6 6 



i=0 j=0 



L-l L-l 



+ E E / / ^( ff l. ff 2) d(7 l (l(J 2^^% 1 ^&-l,fi-l! 



i=0 j=0 



2tt 2tt 



o(— ) + ^2 y y i/j(a 1 ,a2)da 1 da 2 x 
o o 



2tt 2tt 



do\do2 



vO o 



(sin 2 (o-i/2) + sin 2 (a 2 /2)) ; 



where the following relation was used: 



+ o 



2tt 2tt 



n 



2-2A N 



i-0 j-0 o 



do\du2 



'sin 2 {a 1 /2) + sin 2 {a 2 /2)] 



A ' 
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Without loss of generality one may assume k = 1, 1 = 1 in the previous equation. Let us estimate 

2tt 2tt 



47T 
72" 



2 L-l L-l 



< 



EE^-o. 0,0)- f j 



L-l L-l 

EE' 



do\do 2 



{sin 2 (a 1 /2) + sin 2 (a 2 /2)) 



< 



i=0 j=0 



& ft 



(rf(te)/2)+rf(fe)/2))' 
1 



(sin 2 (a 1 /2) + sin 2 (a 2 /2)Y 



do\do 2 



+ 



ft ft 







{sin 2 {ai/2) + sin 2 {a 2 /2)) 



-do\do^ 



+ 



ui + u 2 , 



where E^E^' means summation over 7^ (0,0). 
Let us estimate ui and u 2 . One has: 



L-l L-l 

/ 



=0 i=o 



5, ft 



(sin 2 fe/2) + sin 2 (a 2 /2)) 



( S m2(( f r 1 )/2)+ S m2((a 2 )/2)) A 
do\da 2 



+ 



+ 



EE' 



=0 j=0 



ft ft 



1 



( S m 2 ((&)/2) + S m 2 ((a 2 )/2))' 



(sin 2 (&/2) + sm 2 (&/2)) A 



do\do 2 



nil + «i2- 



The expressions tin and ui 2 can be estimated similarly. Let us estimate u\\\ 

L L 



^EE' 



< 



5^ (^n 2 (fe)/2) + S m 2 ((0)/2)) 



l+A 



< 



L L 



EE' 



< c- 



J2-2X ^2_|_j2)l+A — J^2-2X' 



where c > stands for various estimation constants. Hence 



U\ < 



L 2-2\ ■ 



Let us estimate u 2 : 



U2 



ft ft 







(sin 2 (CTi/2) + sm 2 (cr 2 /2)) 



-da d o 2 



< 
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Si Ci 



< c 



Using polar coordinates, one gets: 



■u 2 < c 



-do\do2- 



l/L 2tt 



L 2-2A • 







Thus: 



U < 



L 2~2\- 

From Lemmas 4.4, 4.1, and Theorem 4.1 it follows that 

2tt 



J 1pl((Tl)d(Tl 



> 







(l + (l))(2vrr+ 1 /^ rg (l) 
2 r r\(rq + l) 1 /?^ - 1 + [R rq (l)]^y ' 



where R rq (t) is a polynomial of degree r, least deviating from zero in L q ([— 1, 1]). 
Theorem 5.4 follows from inequalities (5.4) and (5.5). 



(5.5) 



5.2. Optimal cubature formulas for calculating integrals (1.1). 
Holder class of functions. 

Let x k := 2kir/n, k = 0, 1, . . . , n, A kt = [x k , x k +i,xi,xi +1 ], k, I = 0, 1, . . . , n - 1, x' k = (x k+1 + 
x k )/2, k = 0, 1, . . . ,n - 1, and (si,s 2 ) G A,j, i, j = 0, 1, . . . , n - 1. 
Calculate the integral Kf by the formula: 

^/ = EE M, *i) / / ; 12 , — + (5-6) 



A fc; (^™ 2 (^) + sin 2 (^)) 

Theorem 5.5. Xei \& = H aa (D),0 < a < 1. Among all cubature formulas (5.1) with p\ = P2 = 0, 
formula (5.6), which has the error 

= (2+^i)h 

1 + a Vn/ 

is asymptotically optimal. Here 7 is defined in (5.1'). 

Proof. Using the periodicity of the integrand, we estimate the error of cubature formula (5.6) 
as follows: 



|-Rnn| ^ 



n—1 n—1 



EE 



/(*',, ^)-/(^,*') 



sm' ~ 2 ' + sm z — 7^- 



2gi~-£j 



do\doi 



< 
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+ 



n— 1 n— 1 



fc=0 Z=0 




EE/ /(M.*r 



/(x-,x'))x 



= r\ + r 2 . 



Let us estimate each of the sums ri and r 2 separately. One has: 

/(o-i,o- 2 ) - j{x' k ,x\) 



n < 



E E 

k=i-M l=j-M " £ 



-do\do<i 



+ 



+ 



n— 1 n—l 

EE' 

fc=0 z=o 



= ru + ri2, 



where Y^' means summation over (k,l) such that 
A fci ^ A*, A* = [xi-M,Xi + M+i;xj-M,x j+ M+i],M = [Inn]. 
Furthermore 

do\do2 



r n < — 



£ (sin 2 ^^ + sin 2 ^^) 



A < 



< 



2nM/n 2vr 

r f ^ <ipd0 clogn 



n J J p 

o o 



2 A— 1 j^a-\-2—2X ^ \ ri n 



Estimating ri2, one can assume without loss of generality = (0,0), and get 



w/mr/n 



n—l n— 1 



»"12 < 4 / / (cji((7i) +uj2(cr2))daida 2 y^y^h k i(s 1 ,S2,(Ti,(T2) < 







fe=0 1=0 



n/nir/n 



n—l n—l 







< 



7t\2+a 



<4 I I ((Ji + (J 2 *)d(J 1 d(T2^2^2h kl (s 1 ,S2,<Tl,(T2) < 

k=0 1=0 

n—l n—l 

EE /ifc ^ si '' S2 ' cri ' cr2 ) - 

do\do2 



< 



1 + a \n 

1 + 0(1) 2 /7T 



fc=0 1=0 
2tt 2tt 



1 + a \n 







(sin 2 ^ + sm 2 f )'' 
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Here 



hki(si,s 2 ;ai,a 2 ) = sup h(si, s 2 ; ai, a 2 ) 

(<71,<72)eA fc ( 



Combining the estimates of m and r± 2 , one gets: 

q Vn/ 



1 + 0(1) /7T\<* 

H < -, , - 2(-) 7 



1 + . 



Let us estimate 7*2. To this end we estimate the difference 



r 2 (M)= / / |/(4>^)-/(^>4)| x 



sm 



First, we estimate 



— .So \ A 



) ^sm 



2 — 1 + sm 2 i 



do\do 2 



[sin 



/ I sin 



+ sin 



2 C2— a; j 



' V 



da\da 2 



< 



2+a-2X ' 



n 



The value r 2 (/c, Z) is estimated similarly for |/c — i| < 3 and \l — j\ < 3. 
Let us estimate r 2 (k, I) for other values of k and I. 
One has: 

r 2 (k,l)= I I |/(4,z',)-/(44)|x 



1 



sm 



2£1 _£1 +s - n 2^S2 



sin* — — - + sin* — 2"^- 



< 



n 



2 (Tl— Xj 

|fc-i| 



do\do 2 < 



n 



+ 



n 



(M 2 "- li+ ; i( "- i;) + sin 2 ") 



1+A 



(« n 2 ° i - i+ ; ( "" i) + 



2 "'2-Xj-+S2(s2-^) 



1+A 



< 



< 



|fc-i| 



n 



+ 



K — .71 



n 



< 



\k-i\ 



n 



+ 



n 



n 



(\k-i\ + \l-j\) 



1+a 



— n a+2-2\ 



(\k-i\* + \i-m 



\k - i\* + \l - j\ 2 



< 



1+A 



< 



21 



< 



(\k-i\ 2 + l/-il 2 )( 1+Q )/ 2 



n a+2-2\ - i|2 + |Z - j| 2 ) 

C 1 



< 



< 



n 



a+2-2\ 



(\k-i\* + \l-j 



, 2 Nl/2-a/2+A - 



To estimate r2, one sums up the last expression over k and I. Without loss of generality assume 
= (0,0). Then 



r 2 < 



n 



a+2-2X 



[n/2]+l [ra/2]+l 

16 + 4 E E ' 



fc=o 1=0 



(#2 + /2)A+l/2-a/2 I ' 



where means summation over and / such that k > 3 or I > 3. 

One has: 



[n/2]+l [n/2]+l 

E E ' 



(#2 + /2)A+l/2-a/2 



< 



< A 



fc=0 1=0 
[n/l]+l [n/2]+l [n/2]+l 



fc=3 



fc=3 Z=3 



(#2 + /2)A+l/2-a/2 



< 



< A < 



1, i/ 2A - a > 1; 
logn, i/ 2A — a = 1; 



l-2A+a 



i/ 2A-a<l. 



Hence 



r 2 <^4 < 



n -(a+2-2A) j ^ 2A-a>l; 
ra _1 k>gn, if 2A — a = 1; 
n -1 , i/ 2A - a < 1. 



Thus, if a < 1, then 



r 2 < o(n a ). 
Combining the estimates of n and r 2 , one gets: 

(2 + o(l)) /7T\« 



1 + q \n/ 



Theorem 5.5 follows from the comparison of this inequality with the lower bound of the value 
Cnn[H a ,a(D)], mentioned in the Corollary to Theorem 5.1. I 

Remark. If a = 1, the cubature formula (5.6) is optimal with respect to order. 

The proof of Theorem 5.5 yields also the following result: 

Theorem 5.5'. Let ^ = H aa (D),0 < a < 1. Among all possible cubature formulas (5.1) with 
Pi = P2 = 0, formula 

n—l n—l 



k=0 1=0 



do\do-2 



( S i n 2 (£=21) +sin 2 



"I" Rnm 
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which has the error 

(2 + 0(l)) 7 /7T\« 



Rn 



is asymptotically optimal. 

To apply formula (5.6), one has to calculate the integrals 

W / 7 — — ^ ( "» 

for k, I = 0, 1, . . . , n— 1. Exact values of these integrals for arbitrary values A are apparently unknown. 
Therefore the procedure of numerical calculation of integrals (5.7) should be given for practical appli- 
cation of formula (5.6). 

Let k = i and I = j. Then the integral Jy is replaced by the integral 



7T 7T 

n n 



da 1 da 2 

Pij = / / r , ft > 0, 

n n 

which can be calculated by cubature formulas ( in particular, Gauss quadrature rule ) with arbitrary 
degree of accuracy because the function 



(^sin 2 ^+sin 2 ^-^ +h 

has derivatives up to arbitrary order. The choice of parameter h is discussed in Section 8. 



Let k = i,l j, and 



Let k i, I = j, and 



Let k 7^ i, I ^ j, and 



4tt 2 / . 2 x' l -x' j \ X 



1- ( -:A- 



4vr 2 / . oxi-x'i . 2 x 'l~ x 'i\ A 
The integral if/ is calculated by the formula 

n— 1 n—l 

^/ = EEP*iM^) + ^»(/.P«.^lrf)- (5-8) 

Formula (5.8) is not optimal since it is not exact on constant functions f(x,y) = const. But one 
can estimate the error of this formula: 

n—l n—l 

\R nn (f,Pll,4,y'l))\ <MY J Y.\ I kl-Pll\+ R nnm, 

k=0 1=0 

where M = max \f(x, y)\. 

The values \Im — p* kl \ are easily estimated, and one gets the conclusion of Theorem 5.5'. 
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Classes of smooth functions 

Theorem 5.6. Assume ip £ W r ' r (l). Let ^ = W r,r (l), and calculate the integral Kip by formula 
(5.1) with pi = r — 1, p2 = r — 1, and n\ = n 2 = n. Then the cubature formula 

2tt 2tt 



Kip 



( Pmn{o~l,0-2)daido-2 







(sin 2 {a\ — s\)/2 + sin 2 (a2 — S2)/2) ; 



+ R mn ((p) 



(5.9) 



is asymptotically optimal. 

Before proving Theorem 5.6, let us describe the construction of the spline p> mn . Let x& = 
2/c7r/n, k = 0, 1, . . . , n. Divide the sides of the squares VL = [0, 2ir; 0, 2ir] into n equal parts. Denote by 
Aki the rectangle Am = [2k7r/n, 2(k+l)ir /n; 2l7r/n, 2(/+l)7r/n], k, I = 0, 1, . . . , n — 1. Let (s\, S2) G Aij. 
First we approximate </?(<7i, a 2) as a function of <7 2 , and construct a spline ip n {(T\, cr 2 ) by the following 
rule. Let o\ be an arbitrary fixed number, < o\ < 2tt. On the segments [2kir/n, 2(k + l)ir/n] for 
k 7^ j — 2, . . . , j + 1, one has: 



r-1 



<Pn(0-l,0-2) = ^2 



1=0 



l\ 



((72 - 2kn/n) 1 + BiSVfa, (k + l)/n) 



where 



o(<ri, cr 2 ) := 93(0-1, cr 2 ) - 2^ 7] (02 - 2/cvr/n) . 



The coefficients -Bz are defined by the equation 

(2(fc + i)*/n - o 2 ) r - £ fr44rn V (2vr(A; + " a2r ' _1 = 

= {-l) r Rri (2vr(2fc + l)/2n; ir/n; u 2 ) , 

where R r i(a, h,x) is a polynomial of degree r, least deviating from zero in the norm of the space L 
on the segment [a — h, a + /i]. On the segment [27r(j — 2)/n, 27r(j + 2)/n] the function ip n {o-\,a2) is 
defined by the partial sum of the Taylor series: 

tPn{<ri,o-2) = ip{(Ti,2irj/n) H - (a 2 - j/n) H h 



(r - 1)! 

We define the function (p n n(o~i, 02) by analogy with the function (p n (o~i, cr 2 ). 

Proof of Theorem 5.6. Let (si,.s 2 ) £ ^ij- The error of formula (5.9) we estimate by the 
inequality 

(p((J 1 ,G2) - <Pnn{<7l,&2) 



\Rnn\ ^ ^ ^ ^ ^ 



k=0 1=0 



sin 2 + sin 2 



-da\do2 



+ 



rt-l n-1 

fc=0 /=o 



</?((Tl, CT 2 ) - V3 n n(o"l, 0- 2 ) 



-do\do~2 



-,2 o-j-si 



-,2 CT2-S2 



(5.10) 
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El 
means summation over (k, I) such that i — + — 1 or < ^ n - 1, 

En 
means summation over the other values of (k, I). 

k,l 

Let us estimate each of the sums r\ and r 2 separately. In addition without loss of generality 
assume that ff (ip(<Ti, a 2 ) — <Pnn(&i, cr 2 )) deride ^ 0. Then 



-1 n— 1 



ri ^ e e i^^i'^) - (^1,0-2)1 yy 



do\da 2 



k=o 1=0 



i+l+[(ra-l)/2] j+i+[( n -l)/2] 

r 2 < 4 E E 



J (<j\n 2 £UZ£1 I gin 2 g2n£2 V 
A fei ^ sm 2 + Sm 2 J 



-( r+1 ) ,A<l/2 

-(r+2-2A) iA> l/ 2; 



(5.11) 



(sir 



fc=i+ 2 /=j+2 ( sin 2 ^r 11 + sin 2 , A 



^ 1p((Tl,CT2)d(Tld(J2- 



i+l+[(n-l)/2]j+l+[(n-l)/2] 

- 4 E E 

fc=j+2 i=j+2 



( sin 2 + sin 2 



) (sin 2 ^±pi + sin 2 ^LL 



^ V (o"i, o" 2 )do-i dcr 2 = r 2 i + r 22 , 



(5.12) 



where ^(<7i,<7 2 ) = (p(ai,a 2 ) - ^(^1,^2) 

V> + (0-1,0-2) = < 



■0 (0-1,0-2) = < 



^(0-1,0-2) , if ip(ai,a 2 ) > 

, if V(o-i,o- 2 ) < 0; 

, if^(<7i,<7 2 ) >0; 

-V>(0-l,0- 2 ) , if 1p(<Tl,(T2) < 0. 



It is obvious 



i+l+[(ra-l)/2] j+i+[( n -l)/2] 

4 E E 



(sh 



l + o(l) 



2tt 2tt 



,2 x k -s ± , • 2 ^;-S2 



k=i+ 2 1=3+2 ( sin^ J-fc 2 ° 1 + sin — 



Let us estimate the integral 



- S ^ A ^ 47T 2 



do\da 2 







(sin 2 f + sin 2 f )' 



(5.13) 



^ i){p\,a 2 )do\da 2 ^ yy (^(cti, a 2 ) -<£> n (o"i, a 2 ]}da 1 d(j 2 



+ 



+ 



(0-1,0-2) - ^(0-1,0-2)^0-1^0-2 



= n + «2- 



(5.14) 



25 



Since the expressions i\ and i 2 are estimated similarly, we estimate only i\. One has: 



2vr 

i\ — max 

n s! 



This integral is a continuous function of s±, which attains its maximum at a point s* , and 

%l+l 



2ir 

i\ ^ — 
n 



^ / |<^V,o- 2 )| 

r!n J 1 1 



re, 



r-1 



(r-l-i)! 



2tt /" ^ Bij{xi +1 - x{)r\ w-j-i 
< -T / ^+1 -0-2) - V : ti— (arj+i - ^2 J 



dcr 2 = 



2^+1 



27T f , , 4 /7r\r+2 



From inequalities (5.14) and (5.15) one gets: 



i < 



E " (1) 



and 



(r + 1) 
2 + o(l) /Try 



2?r 2tt 







do\da 2 



(sin 2 f + sin 2 f )' 



One has: 



r 22 = o(n r ). 
Estimate (5.17) follows from the inequalities: 



J J 4> (ai, <J 2 )daida 2 ^ JJ \ifj(ai,a 2 )\daida 2 = 0(n r 2 ) 



and 



i+l+[(n-l)/2] j+l+[(n-l)/2] 

E E 

fc=i+2 i=j+2 



sju 2 - + siu 2 - T ^) A (sin 2 2^±i^£i + S in 2 2i±l=£a^ 



a»TF (fc ~* ) + (/ ~ j) A+ i ^ 
{{k- l y + {i- 3 y) x+1 



n ,A<l/2 
rologn , A = 1/2 
n 2A ,A>l/2. 
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The estimate 



, x , , ^^ 2TT r R rl (l) 

*-<•> 5 (1 + °< i » ( r+ i).(„- 1+ [ fl ; l(1 )| 1 A)^ 



follows from inequalities (5.10), (5.11), (5.16), and (5.17). 

Theorem 5.5 follows from the comparison of the values Cnn[^] and R nn [if]. ■ 

Let us construct cubature formulas for calculating integrals Kf on classes of functions W rr {\). 

These formulas will be less accurate than the ones in Theorem 5.3, but they will be optimal with 

respect to order, and easier to apply. 

First, we investigate the smooth function 

f(Tl,T 2 )dTldT2 



Wuh) = f J 



- - i sin — o 

V 2 



2 Il=ti + sin 2 2i=t2 ) 



A 



assuming f(t±,t 2 ) G W r ' r . Changing the variables t\ = t\ — t, t 2 = t 2 — t, in the last integral, one 
gets: 

f(n +h,r 2 + t 2 )dTid,T 2 



^(tl,t 2 ) = 



Thus, i/j{h,t 2 ) e w r > r . 



(sin 2 ^ + sin 2 a) 



A 



Remark. It is known [9] that Kolmogorov and Babenko widths on the class of functions W r,r (l) 
are equal to 5 n (W r ' r (l)) x d n (W r ' r (l), C) x ^772- Hence the recovery of the function ip(ti,t 2 ) using n 
functionals is not possible with accuracy greater than O (y^rpi^ ■ More precise conclusions are obtained 
in Theorems 5.3 and 5.4- 

Thus, for recovery of a function ip(ti,t 2 ), (ti,t 2 ) £ [0,27r] 2 with the accuracy 0(n~ r / 2 ), it is 
sufficient to calculate the value of the function il>{t\,t 2 ) at the nodes (vk,vi), where Vk = 2kir/N, 
k, I = 0, 1, . . . , N, and iV 2 = n, and to use the local spline ipN{h,t 2 ) of degree r with respect to each 
variable. 

Let us describe the construction of such spline. 

Assume for simplicity that M := N/r is an integer, and cover the domain [0, 27r] 2 with the squares 
Afci = [wk, wi], k,l = 0, 1, . . . , M — 1, here Wk = 2kir/M, k = 0, . . . , M. Approximate the function 
ip(ti,t 2 ) in each domain by the interpolation polynomial ipN{ti, t 2 , A^) constructed on the nodes 
h j = 0, 1, • • • , r, x\ = w k + i = 0, 1, . . . , r. 

Denote the local spline, which is defined by the polynomials ipN(ti,t 2 , A^i), by ipN(ti,t 2 ). 

If the values ip(vk,vi) are calculated by formula (5.9) with the accuracy 0(n" r / 2 ), then 

\mi,t 2 ) - ^ N (ti,t 2 )\\c < 0(n~ r / 2 ). 

Therefore the spline ipN(ti,t 2 ) is optimal with respect to order, and a method for recovery of the 
function ip(ti,t 2 ), which has the error 0(n~ r / 2 ) (in the sup —norm) is constructed. 
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6. Optimal methods for calculating integrals of the form Tf. 
Lower bounds for the functionals ( mn and £jy. 

First we get a lower bound for the error of formula (2.1) with p\ = p2 = and n\ = ri2 = n, on 
Holder classes. 

Theorem 6.1. Let ^ = H aa (D), and calculate the integral Tf by formula (2.1) with n\ = n2 = n 
and p\ = p2 = 0. Then the estimate: 

c m> (1 + 0(1)) // dtldt2 (en 

holds. 

Proof. Let n > be an integer, L = [n/ logn]. Let vu ■= —l+2k/L, k = 0, 1, . . . , L. By (£k,Vl) we 
denote a set which is the union of nodes (xi, tjj), i,j = l,2,...,noi formula (2.1) and the nodes (vi,Vj), 
i, j = 1,2, . . . ,L. Let A kl = [v k ,v k+1 ;vi,v l+1 ], k, I = 0, 1, . . . , L- 1. Let < ip(h,t 2 ) G H aa (D), where 
D = [—1, l] 2 , vanishing at the nodes r)i), k, I = 0, 1, . . . , N, where N = n + L. 

Consider the integral 

l l 

{Tip)(vi,Vj) - 



((ri-^) 2 + (r 2 -^) 2 ) A 

l -l 

'L-lL-1 L-lj-1 i-1 L-l t-1 

EE+EE+EE+EE 

fe=i l=j k=i 1=0 k=0 l=j k=0 1=0 



11 



1p{T 1 ,T2)dT 1 dT2 



((ri-^) 2 + (r 2 -^) 2 ) A 

L-i-lL-j-1 , . 2A 



>- e e d) ((i+1)>+(l+1)J)1 y y «»•»)*>*>+ 

L—i—lj — 1 / \ 2A f f 

+ g gUJ ( (fc + 1)2 + ( z + 1)2^/ / ^i,r 2) ^r 2+ 

A fe+l ,j_i_i 

i-l L-j-1 , L ^2X 1 - - 

+EEUJ ((jfe+1)2+(z+1)2) A y y ^1,^^+ 



fc=0 Z=0 



-fc-lj+i 



^i-fc-lj-i-l 

fc=0 Z=0 v 7 



((fc + l)2 + ( Z + 1 )2)A 



fc=0 Z=0 



28 



k+i,j-l-l 



L 11 ' 'L\ 2X U(i-l-k)U(L-j-l-l) 



+EE 

k=0 1=0 
L-l L-l 

+EE 



((fc + 1)2 + (J + 1)2)A 



i>(jl,T2)dTldT2 + 



L\, 2A - 1 - fc)t7(j - 1 - /) 
((fc + l) 2 + (Z + l) 2 ) A 



y y ip(Ti,T 2 )dTidT 2 . 



&i-k-l,j-l-l 



k=0 1=0 

Here f/(/c) = 1 for jfe > 0, and U(k) = for k < 0. 

Averaging the above inequality over all i and j, i, j = 0, 1, . . . , L — 1, one gets: 

j L-l L-l 



i=0 j=0 



L-l L-l 



> 



L 2-2A 2 2A E E ((jfe + 1)2 + (J + 1)2)A 

xU(L-j-l-l)J y V(n,r 2 )dridr 2 + 



L-l L-l 
i=0 j=0 



L-l L-l 



+ Y,Y, U ( L - i ~ 1 - k ^ U ti - 1 - / / ^(ri,^)drid75+ 

i=o i=o ■'a ^ 

L-l L-l . . 

+ Yl E " 1 " - J - 1 - / / V(Ti,r 2 )dridr 2 + 



1=0 jr=0 

L-l L-l 

:e 

=0 j=0 



^i-fc-lj+i 



L-1L-1 . . 

+ Y, E U ^ - 1 - k ^ U ti - 1 - *) / / ^(n, r 2 )dndr 2 



L-l L-l 



1 1 



> 



2Ao2A 



y y tp(T 1 ,T 2 )dT 1 dT 2 + 



VL-k-2 1 

+y y ip(Ti,T 2 )dT 1 dT 2 + y y ^(^1,^)^1^+ y y ^(n,-^)^^ 



1 



Pk VI 

VL-k-2 Vl-1-2 



Vk -1 



-1 «->! 



-1 -1 



> 



L-l L-l 



1 1 



> 



2Ao2A EE 



L 2-2A 2 2A 



k=0 1=0 

From inequality (6.2) it follows that 



((fc + 1) 2 + (/ + !) 



2\A 



y y i>(Tl,T 2 )dT!dT 2 . (6.2) 



-1 -1 



L-l L-l 



1 1 



Cr, 



>[Haa{D)] > (1 + 0(1)) E E ( fe 2 + /2)A / / ^T^dr^ 

k=l 1=1 



2 2A 4 



1 1 



1 1 



S S w^w S S ^ T2)dTldT2 - 

-1-1 -1-1 



(6.3) 
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From Theorem 4.2 and Lemma 4.4 it follows that the inequality 

1 l 

/ / 4>{Tl,T 2 )dT 1 dT 2 > -4 \ (6-4) 

J J l + an a 

-i -l 

is valid for an arbitrary vector of the weights and the nodes (X,Y,P) on the class H aa (D). 
Theorem 6.1 follows from inequalities (6.3) and (6.4). ■ 

Theorem 6.2. Let ^ = C^l), and calculate the integral Tf by formula (2.1) with p\ = P2 = 0. 
If n\ = 7i2 = n, then 

l l 

2K r f f dsids2 



c™ffl>a+°a)>^// 



-i -i 

where K r is the Favard constant. 
Proof. Let 

i>{si,s 2 ) = tpi(si) + ^2(52), 

where < ipi(s) G W r (l), vanishes at the nodes Xk, k = 1,2, ... ,n, and < ip2(s) £ W r (l) vanishes 
at the nodes y^, k = 1,2, ... ,n. 

For arbitrary nodes x^, k = 1, 2, . . . , n, one has (see [11]): 

1 

/2K r 
A(s)ds > - — i = 1,2. 
(7rn) r 

-1 

Thus the inequality 

J J i>(s 1 ,S 2 )ds 1 dS2 > 

-1 -1 

holds for arbitrary nodes {x\, . . . , x n ) and (yi, . . . ,y n ). 

Theorem 6.2 follows from this estimate and inequality (6.3). ■ 

Theorem 6.3. Let ^ = Wp T {l), r = 1, 2, . . . , 1 < p < 00, and calculate the integral Tf by 
formula (2.1) with p\ = p2 = r — 1 and n\ = ri2 = n. Then the estimate 

1 1 

, r 1 , , „ 2 1 /iR ra (l) f f d Sl ds 2 , s 

-1 -1 

holds, where R rq {t) is a polynomial of degree r, least deviating from zero in L q ([—1, 1]). 

Proof. Let L = [n/logn]. Consider the nodes (vk,vi), Vk = k,l = 0,1,... ,L — 1. By 
(Cii^j)) = 0, 1, . . . , AT — 1, iV = n + L denote the union of the nodes (xk,yi) and (£«,£•,). Let 
^(si,s 2 ) = ^i( s i) + ^2(^2), where < Y>i(s) £ Wp"(l) vanishes with its derivatives up to order r — 1 
at the nodes z = 0, 1, . . . , N — 1, and < "02 (s) £ Wp (1) vanishes with its derivatives up to order 

Vi+l 

r — 1 at the nodes r/j, j = 0, 1, . . . , N — 1. Assume that f ip\(s)ds > 0, i = 0, 1, . . . , N — 1, and 
/ ^(s)ds>0, j = 0, 1, . . . , AT — 1. 
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Using the argument similar to the one in the proof of Theorem 6.1, one gets: 

j L-1L-1 
i=0 j=0 

L-1L-1 1 1 



i^X E E ((fc + 1)2 + (l + 1)2) A / / ^(n, r 2 )dndr 2 



fc=0 / _ , _ , 

11 11 
1 + o(l) /" /* dtidt 



y y ^{T 1 ,T2)dT 1 dT 2 . (6.6) 



2 2A 4 y y (tf + ti) A 

-l-i -l-i 
Prom Theorem 4.1 and lemma 4.4 it follows that the inequality 

//*<™>*>*> > ,1 + °m) rKrq+iy ilZ^it,mv>°r <6 ' 7) 

is valid for arbitrary weights and the nodes (X, Y, P) on the class H aa (D). 
Theorem 6.3 follows from inequalities (6.6)- (6.7). ■ 
Cubature formulas. 

Let us construct a cubature formula for calculating the integral Tf on the Holder class 
Haa(D). Let x k '■= —1 + 2k/n, k = 0,1, ... ,n, x' k = (x k +i + x k )/2, k = 0, 1, . . . ,n — 1, and 
Afci = [x k ,x k+ i;xi, x i+1 ], k, I = 0, 1, . . . , n - 1. 

Calculate the integral Tf by the formula 

n— In— 1 „ „ , , 

r/= 2S /(4, l!) / / (fr + <6 - 8) 

Consider another cubature formula for calculating the integral Tf. 

Let (ii, t 2 ) € Ay. By A* denote the union of the square Ajj and of those squares A^ which have 
common points with the Ay. Consider the formula 

T/ = /( ^' ^ / / ((ri-ti) 2 + (r 2 -t 2 ) 2 ) A + 
A* 

n— In— 1 „ „ , , 

+ EE'/fex!)// ((Ti _ ti) y;_ t2)2)> + ^(/). (6.9, 

where ^ ^ ' means summation over the squares which do not belong to A* . 

Theorem 6.4. Among all cubature formulas (2.1) with p\ = p 2 = and n\ = n 2 = n, formula 
(6.8), with the error estimate (6.15), is optimal with respect to order. 

Remark. Similar statement holds for formula (6.9). 

Proof of Theorem 6.4. Let us estimate errors of formulas (6.8) and (6.9). 
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The error of formula (6.8) can be estimated as follows: 

n ~ ln -\ " - _|/(r 1 ,r 2 )-/(4^ 

A---U /— (> 



iu/)i<EE' 



(( Tl _t 1 )2 + (r2 _ t2) 2 ) A 



dT\dT2 + 



n—1 n—1 

+EE" 

fc=0 z=o 



(( Tl _ tl )2 + (r2 _ t2) 2 ) A 



dT\d,T2 = r\ + r 2 , 



(6.10) 



where YJ YJ' means summation over fc and / such that the squares A&; belong to A*, and YJ YJ" means 
summation over the other squares. 
Let us estimate r\ and r 2 : 



n - n«J J ((n - ti) 2 + (r 2 - i 2 ) 2 ) A " n 2 ^+« 
A* 

n—1 n—1 



= o(n- a ); 



(6.11) 



(6.12) 



fc=0 «=o 



Here h(Aki) denotes the maximum value of the function ((n — ii) 2 + (r 2 — t 2 ) 2 ) A in the square Aki- 
One has: 

1 



L((ri-ti) 2 + (r 2 -t 2 ) 2 ) 



2U 



kij 



dr\dT2 



A fe ; 
< 



L ((n - hf + (r 2 - t 2 ) 2 ) A ((x k - + (x t - t 2 ) 2 ) x \ 
2\(x k -h + qi(n - x k ))(ji - x k )) 



II 



+ 



((x k ~h+ ft (n - x k )f + (r 2 - t,) 2 )^ 1 
2A(x; - t 2 + 9 2 (r 2 - x/))(r 2 - x/)) 



dr\dT 2 



dT\dr 2 



+ 



< 



((zfc - *i) + (xi -t 2 + q 2 (r 2 - x t )) 



2\X 



dT\dr 2 



< 



<-II 



2A(n - x k ) 

((x fc - h + g(n - x fc )) 2 + (r 2 - t 2 ) 2 ) A +V 2 



dr\dT 2 



+ 



2A(r 2 - a;)) 

((n - ti) 2 + (s, - * 2 + g 2 (r 2 - x ; )) 2 ) A+1 / 2 



dT\dr 2 



+ 



< 



< 



2 4 A n 2A+1 



2 4 A 



1 



n 3 ( fc 2 + /2)A+l/2 n 2-2A (fc2+/ 2)A+l/2' 

where it was assumed that k > i + 1, and / > j + 1. Estimates for the other combinations of k and I 
are similar. Thus: 

n—1 n—1 „ „ , , 

OTi GET 



1«T 2 



r2 "(l+«)^S"i 7 ((TL-tl) 2 + (71,-t2) 2 ) A 
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1 2 3 A n - ln - 1 i 
+ (l + a) n 2 ' 2X + a (fc2 + /2)A+i/2- ( 6 - 13 ) 



Let us estimate the last term in the above inequality. 
One has: 

n-ln-l 1 n/2 ra/2 

VV" - < V V * - < 

to h ^ + ' 2 ) A+1/2 - J( m J&n ^ + l2 ) x+1/2 ~ 

1, A>l/2 

c < logn, A = 1/2 (6.14) 
n l ~ 2X , A<l/2 
where ^ YT means summation over k and /, (k,l) / (0,0). 

In deriving (6.14) we have used the known result ([17], Theorem 56) which says that a number 
of points with integer- value coordinates, situated in the circle x 2 + y 2 = r 2 , is equal to nr 2 + 0(r). 
From inequalities (6.13) and (6.14) it follows that 

(l + o(l)) 1 ^^\, f f d n dr 2 



z — fU^ n« ^ 

This and (6.11) yield: 



+ na toU J J ((n-h) 2 + (r 2 -t 2 ) 



l l 

dT\dr 2 



2\A • 



Rnn[H aa (D)}< ( ]^\ SUp / / 

(1 + a)n a {tl ,t 2 )eDJ J 



((ri-ti) 2 + (r 2 -t 2 )2)A 

l -l 



< 



< ; +oW (6.15) 

- (1 + a)n a J J (t 2 + r|) A V ' 

-l -l 

Theorem 6.4 follows from a comparison the estimates of C[H a ^ a (D)} and R nn [H ctja (D)]. ■ 
Let us construct optimal with respect to order cubature formula for calculating integrals Tf 
on the classes W rr . In the derivation of formula (5.9) the local spline (p n (ti,t 2 ), approximating the 
function <p(t±,t 2 ) in the domain [0, 2tt; 0, 2tt], was constructed. A spline f n n{ti,t 2 ), approximating 
the function f(t±,t 2 ) in the domain [—1,1] x [—1,1], can be constructed analogously. Calculate the 
integral Tf by the formula 

-i -i 

Theorem 6.5. Let ^ = W r ' r (l), r = 1, 2, . . . , and calculate the integral Tf by formula (2.1) with 
pi = P2 = r — I; and n\ = n 2 = n. Then cubature formula (6.16), which has the error 



l i 

J? ATA ft _L_ fU\ 2 ^rl(l) f f d Tl dT 2 

-1-1 
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is optimal with respect to order. Here R rq {t) is a polynomial of degree r, least deviating from zero in 
As in the proof of the Theorem 5.6 one gets the following estimate 

*-<•> £ (1 + * \ + m n^Lwn 1 1 

-1-1 

Comparing this estimate with the estimate of Cnn[W r ' r (l)] from Theorem 6.3 one finishes the 
proof. ■ 

7. Calculation of weakly singular integrals on piecewise continuous surfaces. 

In Sections 5 and 6 asymptotically optimal methods for calculating weakly singular integrals 
defined on the squares [0,2tt] 2 or [-1,1] 2 were constructed. 

It is of interest to study optimal methods for calculating weakly singular integrals on piecewise- 
Lyapunov surfaces. 

Consider the integral 

Jf = [f /(Tl ' r2 ' T3)d5 ~ x , h,t 2 ,t 3 e G, (7.1) 

G ((Ti-ti) 2 + (r 2 -t 2 ) 2 + (T3-t 3 ) 2 ) 

where G is a Lyapunov surface of class L s (B,a). 

We show that the results derived in Sections 5 and 6 can be partially generalized to the integrals 
(7.1). 

Calculate integrals (7.1) by the formula: 

n p 

i/ = EE Pkvf {v \M k ) +Rn(f,G,M k ,p kv ,t), (7.2) 

k=l \ v \=0 

where t = {h,t 2 ,t 3 ), v = (vi,v 2 ,V3,), \v\ = v\ + v 2 + v 3 , f^(h,t 2 ,t 3 ) = q^^Iq^ ■ 
The error of formula (7.2) is: 

Rn(f,G,M k ,p kv ) = swp\R n {f,G,M k ,p kv ,t)\. 
tec 

Assume / € Vl/i, and G € VP 2- Then the error of formula (7.2) on the classes Vti and Vt^ is: 

R n (* 1 ,V 2 ,M k ,P kv )= sup R n (f,G,M k ,p kv ). 

fev u Ge*2 



Let 



Cn[*i,*2] := „mf Rn^V^MfrPto). 

M k ,p kv 



A cubature formula with nodes and weights p* kv is called optimal, asymptotically optimal, 
optimal with respect to order on the class of functions vPi and surfaces ^> 2 , if 
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respectively. 

Let = H a (l), < a ^ 1, and $>2 = L\{B,(5) < (3 ^ 1. Let us construct an optimal with 
respect to order method for calculating integrals (7.1) on the classes of functions ^1 and surfaces \l/ 2 . 
Let S(G) be a "square" of the surface G. Divide the surface G into n parts gk, k = 1, 2, . . . , n, so that 
a "square" of each of the domains g\. has the area of order \S(G)\/n, where \S(G) | is the area of S(G). 
We take a point in each of domains gk at the center of the domain gk- 

Calculate integral (7.1) by the formula 

dS 

f{M k ) J J 



Jf = J2f(M k ) 



n ((n - ti) 2 + (r 2 - t 2 y + (T3 - t 3 y) 

Theorem 7.1. Formula (7.3), has the error 



+ Rn(f,G). 



(7.3) 



and is optimal with respect to order on the classes *$>i = H a , < a ^ 1, and ^2 = L\{B,(3), < (3 ^ 
1, among all formulas (7.2) with p = 0. 

Proof. Assume for simplicity that the surface G is given by the equation z = ip(x,y), (x,y) G 
Gq, tp(x,y) ^ 0. Let (p x (x,y) := p, ip y (x,y) := q. Write the integral Jf as 



Jf 



f(n, T2, <p(ti,t 2 )) v / 1+p 2 {ti,t 2 ) + q 2 (Tl,T 2 )dT 1 dT 2 



Go 



(ri - ti) 2 + (r 2 - t 2 ) 2 + Mri,r 2 ) - ip(h,t 2 )) 2 



x • 



(7.4) 



The function /(n, r 2 , v?(ri, r 2 )) belongs to the Holder class H a over Go, and the function 

y / l+p 2 +9 2 



is positive. 

(ri-ti) 2 + (T2-t2) 2 + (v(Ti,r 2 )- V .(t 1 ,t 2 )) 2 

Let Mk = {m\, m\,m\) be the nodes of cubature formula (7.2). Let iP(t) := (d(r, {M^}))", where 
<i(r, {M/J) is the distance between the point r and the set of the nodes {M^}, where the distance 
is measured along the geodesies of the surface G. This distance satisfies the Holder condition H a {\). 
Hence the function ip*(T\,T2) = 4>(ti,T2,ip(ti,T2)) belongs to the Holder class H a (A) and vanishes at 
the nodes {m!{, m 2 ), k = 1, 2, . . . , n. Thus, 



Cn(*l,*2) > 



S(G ) 



x mm 

t 



G 



Go Go 



V ; (n,T 2 ,y(ri,r2))v / l +p 2 + q 2 dT 1 dT 2 dt 1 dt 2 
(ri - h) 2 + (r 2 - t 2 ) 2 + (y(ri,r 2 ) - <p(ti,* 2 )) 2 



^ S(G ) / / ^( Tl ' r 2'^( r l' r 2)) (ir l d ' r 2 X 
Go 



y 7 ! + p 2 + Q 2 



(ri - ti) 2 + (r 2 - t 2 ) 2 + Kri,r 2 ) - <^i,i 2 )) 2 



- 1 -dridT2 ^ 



^ — mm 



// ds 

na/2 1 JJ (r(t,r)) X ' 
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where S(Go) is the "square" of the surface Gq. 

Therefore the error of formula (7.3) is estimated by the inequality R n ^ 



n' 



A 



Theorem 7.1 is proved. ■ 

Remark 1. The method of decomposition of the domain G into smaller parts gk, fc = l,2,...,n, 
described below, is optimal with respect to order for classes of functions = H a , < a ^ 1, and of 
surfaces * 2 = L (B, /3), < (3 < 1 for a < (3. 

Remark 2. From formula (7.4) it follows that if the function f £ W r,r (l) and the surface G £ 
L s (B,a), then the function /(ri, t 2 , </>(ti, t 2 )) £ W V ' V {A), where v = min(r, s). Therefore, repeating 
the above arguments, one proves that the accuracy of calculation of integral (7.4) by cubature formulas 
using n values of integrand function does not exceed 0{n~ v / 2 ). 

From this remark it follows that if the surface G consists of several parts, for example of surfaces 
G\ and G 2 having common edge L, then it is necessary to calculate the integrals for the surface G\ 
and the surface G 2 separately. If the surface G is divided into smaller parts gk, k = 1,2, . . . , n, the 
domains gk, the curve L passes inside of these domains, should be associated with the class of surfaces 
L (B, 1). In these domains the accuracy of calculation of the integral does not exceed than 0(n^ X ), 
where nk is the number of nodes of the cubature formula used in the domain g^. 

For this reason the cusps and the nodes, in which three or more domains Gk, which are parts of the 
domain G touch each other, must belong to the boundaries of the covering domains gk, k = 1, 2, . . . , n. 

The universal code for computing the capacitances, described in Section 9, is based on optimal 
with respect to order cubature formulas for calculating integrals on the classes of functions H a , < 
a ^ 1 and of surfaces L (B,/3), B = const, a ^ (3, (3^1. 

The algorithm constructed in Section 9 is optimal on this class of surfaces and does not require 
special treatment of edges and conical points of the surface. 

When one studies cubature formulas on the classes W r ' r (A), r > 1, and L S (B,(3), s ^ 1, 0^/3^ 
1, one has to develop a method to compute accurately the integrals in a neighborhood of the above 
singular points of the surface. 

8. Calculation of weights of cubature formulas. 

In calculating weakly singular integrals by cubature formulas (6.8) it is necessary to calculate 
integrals of the form of 



for different values (ti,t 2 ) £ [— 1, l] 2 . 

Let (ti,i 2 ) £ Ay. Let us consider two possibilities: 1) the square A/y and the square Ajj have 
nonempty intersection; 2) the square Aki is does not have common points with the square A^-. 
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First consider the second case, when the function 



l 



((ri-ti) 2 + (r 2 -t 2 )2) A! 

is smooth. Here (ti,T2) G A^, and (ii,t 2 ) G A^. 
In this case one has: 



d r (f(Tl,T 2 ) 



< 



, 2 2r 



(( ri _ tl )2 + (r2 _ t2)2) A+r/2 



and, if the squares A^ and Ajj do not have common points, one gets: 



9 r V?(ri,r 2 ) 



< 



2 r r\n 2X+r 
~2~ x • 



Similar estimates holds for partial derivative with respect to r 2 . 
Calculate the integral Jki{ti,t2) by the Gauss cubature formula: 



Jki(h, h) — j 



((ri-ti) 2 + (r 2 -t 2 ) 



2^ A 



dndT2 + Rmm(^kl), 



where P m m = Pm^m > Pm (* = 1) 2) is the projection operator onto the set of interpolation polynomials 
of degree m with nodes at the zeros of the Legendre polynomial, which maps the segment [—1,1] onto 
the segment [xk,Xk+i] for i = 1, and onto the segment [x/,x/ + i] for i = 2. 



An integer m is chosen so that \R r , 



< n 



-2-a 



for cubature formulas on the Holder class H n 



and |-R mm | < n" r ~ Q for cubature formulas on the class W rr . 

This requirement is made because the error of calculation of the coefficients Jfc/(ti,t 2 ) must not 
exceed the error of formula (6.5). 

Using r derivatives of the integrand in the error R mm (Aki), one gets: 



\Rmm{^kl)\ < 



B r 2 r r\ (2 



,r-l 



2-2A 



m' x \n / 

where B r is the constant appearing in Jackson's theorems. It is known that the constants B r are 
bounded by a constant, denoted b, uniformly with respect to r. In the case of periodic functions 6 = 1 
([18]), and in the general case b is apparently unknown. 

If r = 2 and m = B r 2 r r\n 2X , then one gets the error estimate given for cubature formula (6.5). 

Now, consider a method for calculating the integrals Jki(ti, £ 2 ) when the square A^; has nonempty 
intersection with the square A^-. For defmiteness we consider the calculation of the integral Jij(t\,t2) 
by the formula: 

1 



Ji 



i(h,t 2 ) = J 



{{Ti-hY + i^-^Yf + h 



dTid,T 2 + Rmm(Aij), 



where h = const > will be specified below. 
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One has: 

\Rmm(\j)\ <h f -p** ^ + 

l ((T! - hf + (r 2 - ttf t (((ri - ti) 2 + (r 2 - t 2 ) 2 ) A + h) 

I 



A, 



((ri-hy + (T2-t 2 y) x + h 



d,T\d,T2 = n + r 2 , 



where D mm = I — P mm , and / is an identity operator, and 



<h f dr\dri 

1 - / fuiy , / „ \(l+A)/2 - 

a4 ((ti - ti) 2 + (r 2 - t 2 ) 2 ) (A+1)/2 (((n - t x f + (r 2 - t 2 ) 2 ) A + h) 

<^ (1 - A)/2 f-V _2A - (8-1) 



1 - A \n J 

((ri-ti) 2 + (r 2 -t 2 )2) A +fe 



The function — — — -r — is infinitely smooth. Using bounds for its first derivatives for 



A > 1/2, one gets: 



r 2 < ^3^. 8.2 



Prom inequality (8.1) it follows that for getting accuracy 0{n~ l ~ a ) one has to have h = 
n -2(2A+a)/(i-A) and from inequality (8.2) it follows that one has to have m = max([n( 8A+4a )/( 1 " A ) +Q " 3 ], 1). 

9. Iterative methods for calculating electrical capacitancies of conductors of arbitrary 
shapes. 

Numerical methods for solving electrostatic problems, in particular, calculating capacitancies of 
conductors of arbitrary shapes, are of practical interest in many applications. Electrostatic problems 
solvable in closed form are collected in [19,20,21]. Some of the problems were solved in closed form 
using integral equations, Wiener-Hopf and singular integral equations [22]. Electrostatic problems 
for a finite circular hollow cylinder (tube) were studied in [23] by numerical methods. In [24] the 
variational methods of Ritz and Trefftz are discussed. Galerkin's and other projection methods are 
studied in [25]. In practice these methods are time-consuming and variational methods in three- 
dimensional static problems probably have some advantages over the grid method. There exists a vast 
literature on calculation of the capacitances of perfect conductors [20,26]. In [20] there is a reference 
section which gives the capacitance of the conductors of certain shapes (more than 800 shapes are 
considered in [20]). In [26] and [27] a systematic exposition of variational methods for estimation of 
the capacitances is given. In [28] there are some programs for calculating the two-dimensional static 
fields using integral equations method. In monograph [7] iterative methods for solving interior and 
exterior boundary value problems in electrostatics are proposed and mathematically justified. Upper 
and lower estimates for some functionals of electrostatic fields are obtained in [7]. Such functionals 
are the capacitances of perfect conductors and the polarizability tensors of bodies of arbitrary shape. 



38 



These bodies are described by their dielectric permittivity, magnetic permeability and conductivity. 
They can be homogeneous or flaky. The main point is: these bodies have arbitrary geometrical shapes. 

The methods, developed in [7], yield analytical formulas for calculation of the capacitances and 
polarizability tensors of bodies of arbitrary shapes with any given accuracy. Error estimates for 
these formulas are obtained in [7] . We give here the formulas for calculating the capacitances of the 
conductors of arbitrary shapes [7]: 



C(«) = 4Tre S 2 



(2 



% I [—[...[ 4>{t,t l )...^(tn-ut n )dt l ...dt n 

tt) h Jr r st Jr Jr 



> 



71 times 

where S is the surface area of the surface T of the conductor, eq is the dielectric constant of the 
medium, r st := \s - t\, and tp(t,s) := 



C (0) 

It is proved in [7] that 



4tt£ S 2 



J 



<C, J 



dsdt 
r s t 



S = measr. 



c - c {n) 



< Aq n , < q < 1, 



where A and q are constants which depend only on the geometry of T. 

We use these formulas are used to construct the computer code for calculating the capacitances 
of the conductors of arbitrary shapes. 

It is proved in [7], that 

C (n) = 4Tre S 2 J J r^ 1 S n (t)dtds 



(9.1) 



r 



where 5 n is defined by the iterative process: 

<Wi = ~A5 n 



, 5 = 1, J S n dt = S, 
r 



(9.2) 



and A is defined by the formula: 



A5 



j m . 



d 1 



-dt, 



8N S 2-nr st 
r 

where N s is the outer unit normal to T at the point s. 

To use iterative process (9.2), one has to calculate the weakly singular integral 

d 1 



dN s r st 



-dt. 



(9.3) 



Let us describe the construction of a cubature formula for calculating integral (9.3), assuming for 
simplicity that the domain G, bounded by the surface T, is convex. This asumption can be removed. 
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Let 8 be the inscribed in the conductor sphere of maximal radius r* , centered at the origin. 
Introduce the spherical coordinates system {r,<p,6), and the set of the nodes (r* ,(f) k ,9{), where (ft k = 
2kir/n, k = 0, 1, . . . , n, B\ = nl/m, I = 0,1, ... ,m. Assume that m is even, and cover the sphere § 
with the spherical triangles A k , k = 1, 2, . . . , N, N = 2n(m — 1). 

Let us describe the construction of the spherical triangles. For < < ir/m the triangles 
A k ,k = 1,2,. . . ,n have vertices (r*,0,0), (r*, (f) k -i, 0i), (r*,(f> k ,9i), k = 1,2, ... ,n. 

For 6i < 6 < 6i +1 , I = 1,2, ...,m/2 - 1, the triangles A k ,k = n + 2n(l - 1) + j, 1 < 
j < 2n are constructed as follows. The rectangle [0, 2tt; 6i, 0i+i] is covered with the squares 
Afci = [4> k , 0fc+ii &h ^i+i], k = 0, l,...,n — 1. Each of the squares A^; is divived into two equal 
triangles A^ and A^, k = 0, 1, . . . , n — 1, / = 1,2, ...,m/2 — 1. The spherical triangles Aj^ and 
A^, k = 0, 1, . . . , n — 1, I = 1,2, ... , m/2 — 1, are images of triangles A^ and A^ on the sphere S 

As a result of these constructions the sphere S is covered with triangles A^, k = 1, 2, . . . , N. 

We draw the straight lines through the origin and vertices of the triangle A&, k = 1,2,... ,N. The 
points of intersection of these lines with the surface T are vertices of the triangle A k , k = 1, 2, . . . , N. 
As a result of these constructions the surface F is approximated by the surface T n consisting of triangle 
A k , k = 1,2, . . . , N, and integral (9.3) is approximated by the integral 



We fix each triangle A k , k = 1,2, ... ,N, and associate with it a point T k G A k , k = 1,2, ... ,N, 
equidistant from the vertices of the triangle A k , k = 1,2,...,N. We calculate integral (9.4) at the 
points T k , k = 1,2, ... ,N, by the cubature formulas constructed in paragraphs 5-7 for the Holder 
classes. After calculating the values U (r^), k = 1, 2, . . . , N by these cubature the integral 



is calculated, where U(t) = U{r k ) for t 6 A^, k = 1,2, ... ,N, Sn is area of the surface Tn, C^> is 



10. Numerical examples. 

In this section the numerical results are given. As an example we calculated the capacitances of 

various ellipsoids, because for ellipsoids one knows ([19]) the analytical formula for the capacitance, 

which makes it possible to evaluate the accuracy of the numerical results. Consider the ellipsoid: 

x 2 y 2 z 2 

h - — I = 1. 

a 2 b 2 c 2 

It is known [19,20] that the exact value of the capacitance of ellipsoid with a = b is: 




(9.4) 




approximation to the value of . The successive iterations are calculated analogously. 



C = 



47r£oV(a 2 - c 2 ) 
arccos(c/a) 
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Let a = b = 1, and £q = 1. We have calculate the capacitance C for different values of the 
semiaxis c. The results of the calculations are given in Table 1. 

It is known ([7], p. 43), that the capacitance of a metallic disc of radius a is C = 8aeo, and one can 
see from Table 1, that asymptotically, as c — > 0, this formula can be used practically for the ellipsoids 
with c < 0.001 with the error about 0.005. 
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Table 1 



c 


n 


m 


N 


Exact value 


Error 


Relative error 


Calculation time 


0.9 


40 


30 


2320 


12.144630 


-0.221200 


0.018212 


25 sec 


0.5 


40 


30 


2320 


10.392304 


-0.222042 


0.021366 


25 sec 


0.1 


40 


30 


2320 


8.5020638 


-0.301189 


0.035425 


25 sec 


0.01 


40 


30 


2320 


8.050854 


0.072132 


0.008959 


25 sec 


0.001 


40 


30 


2320 


8.005092 


-0.821528 


0.106374 


25 sec 


0.0001 


40 


30 


2320 


8.000509 


-1.068178 


0.133513 


25 sec 


0.9 


50 


40 


3900 


12.144630 


-0.180510 


0.014801 


1 min 15 sec 


0.5 


50 


40 


3900 


10.392304 


-0.185642 


0.017860 


1 min 15 sec 


0.1 


50 


40 


3900 


8.5020638 


-0.288628 


0.033947 


1 min 15 sec 


0.01 


50 


40 


3900 


8.050854 


-0.372047 


0.046212 


1 min 15 sec 


0.001 


50 


40 


3900 


8.005092 


-0.586733 


0.073295 


1 min 15 sec 


0.0001 


50 


40 


3900 


8.000509 


-0.933288 


0.116653 


1 min 15 sec 


0.9 


60 


50 


5880 


12.144630 


-0.152009 


0.012516 


4 min 


0.5 


60 


50 


5880 


10.392304 


-0.160023 


0.015391 


4 min 


0.1 


60 


50 


5880 


8.5020638 


-0.283364 


0.033328 


4 min 


0.01 


60 


50 


5880 


8.050854 


0.532250 


0.061110 


4 min 


0.001 


60 


50 


5880 


8.005092 


-0.391755 


0.048939 


4 min 


0.0001 


60 


50 


5880 


8.000509 


-0.880394 


0.110042 


4 min 
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